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Abstract 

Webs are sets of Feynman diagrams that contribute to the exponents of scattering ampli¬ 
tudes, in the kinematic limit in which emitted radiation is soft. As such, they have a number 
of phenomenological and formal applications, and offer tantalising glimpses into the all-order 
structure of perturbative quantum held theory. This article is based on a series of lectures given 
to graduate students, and aims to provide a pedagogical introduction to webs. Topics covered 
include exponentiation in (non-)abelian gauge theories, the web mixing matrix formalism for 
non-abelian gauge theories, and recent progress on the calculation of web diagrams. Problems 
are included throughout the text, to aid understanding. 


1 Introduction 

One of the purposes of quantum field theory is the calculation of scattering amplitudes, which are 
related to the probability for a given particle collision to occur. Typically, amplitudes are calculated 
in perturbation theory, and it is well known that such calculations are beset by infinities. Some 
of these infinities (ultraviolet divergences) can be absorbed into redefinitions of the Lagrangian 
parameters defining the theory, a procedure known as renormalisation. In these lectures, we will 
be concerned with infrared (IR) divergences, which correspond to long distance physics. 

That long distance physics is associated with divergences in perturbation theory can be inferred by 
examining figure This depicts a 2—7-2 scattering process, in which one of the outgoing particles 
emits a gluon. This emission may be real, in which case the gluon itself ends up as a final state 
particle. Alternatively it may be virtual, in which case it joins up again with the outgoing particle. 
These two possibilities are shown in figure]^ a) and (b) respectively. The Feynman rules for emission 
of a gluon tell us to integrate over all possible positions for the end of the gluon. That is, there is 
an integral 



in d space-time dimensions. Because one integrates out to infinite distance, there will always be 
some value of d above which this integral diverges, and in fact this is found to be d = 4. Hence, 
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(a) (b) 

Figure 1: Emission of a gluon of momentum k from an outgoing particle in a scattering amplitude: 
(a) real emission; (b) virtual emission. 


singularities associated with long distance physics are a crucial ingredient in practical applications 
of field theory. 


Note that a gluon that extends out to infinite distance has an infinite Compton wavelength. Thus, 
by the uncertainty principle, such a gluon has a vanishing momentum or energy —?• 0. This is the 
origin of the phrase “infrared divergences” for long distance singularities - if we were calculating 
Feynman diagrams such as those in figure in momentum space, the singularities would show up 
after integrating over the emitted gluon momentum k^, as contributions from the region where this 
tends to zero in all components. Such gluons are referred to as soft in the literature, as opposed to 
hard particles which have nonzero momenta. 


From hgure we also see that if a gluon has an infinite Compton wavelength, it is difficult 
to distinguish between the case of real or virtual emission, as a virtual gluon will only rejoin an 
outgoing particle at inhnity. Thus, both real and virtual diagrams will contain infrared singularities 
(after integrating over the gluon position in the real case), and we expect that the real and virtual 
contributions for a given soft gluon will be closely related to each other. This is indeed the case, 
and in fact IR divergences cancel upon combining real and virtual graphs. Furthermore, this is 
guaranteed, for suitably defined quantities, to all orders in perturbation theory [^. However, large 
residual contributions to observable quantities remain, which typically take the form of terms such 


as 


ag- 


, log’" E 
E 


( 1 ) 


involving logarithms of the total energy E carried by soft gluons. Such terms occur at all orders 
0{a'g), and are highly divergent as E —)• 0 such that perturbation theory breaks down in this 
regime. Close to this limit, the logarithmic terms are not formally singular, but can be numeri¬ 
cally large. Then perturbation theory is no longer sufficiently convergent for a truncation at e.g. 
next-to-leading order to be sensible. One must then sum up the terms of eq. ([^ to all orders in the 
coupling constant, a procedure known as resummation. This is by now a highly developed held - 
3.g. [2-12 for several different approaches, and refs. |^,14 for recent review articles. Implicit 


see e. 
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in the above discussion is that resummation is intimately linked with the structure of infrared sin¬ 
gularities to all orders in perturbation theory. 


Infrared singularities are also of interest in themselves. In particular, there are a number of recent 
conjectures regarding IR divergences in a variety of theories. One example is the so-called dipole for¬ 
mula in QCD 15 -19 , an all-order ansatz for the IR singularities of fixed-angle scattering amplitudes 
with massless particles. There is also interesting recent work on singularity structures in super- 


symmetric gauge theories, and their relationship to similar amplitudes in (super-)gravity (20-25 
The IR singularity structure of quantum gravity itself has been studied in , and more recently 
The relation between QCD and gravity in the soft limit has been explored in 


in 27-30 


used to provide evidence for the double copy conjecture relating both theories 20,21 


31 , and 


As a result of these investigations, much is already known about the structure of infrared diver¬ 
gences. Crucially, it is known that scattering amplitudes factorise into the following schematic 
form (^-38 : 

^ J 
Ji' 


A^n s Yl 


( 2 ) 


2=1 


Here the left-hand side represents an amplitude with a fixed number L of outgoing hard particles, 
where it is assumed that ultraviolet renormalisation has already taken place. The first factor on 
the right-hand side is the hard function, and is finite in four dimensions (i.e. contains no IR sin¬ 
gularities). We can think of this as a subamplitude for the outgoing hard particles, where no soft 
gluons are emitted. The second factor S is the soft function, and collects all singularities associated 
with gluon emissions whose 4-momenta each tend to zero. The final factors, one for each outgoing 
particle, collect collinear singularities, i.e. divergences associated with gluons emitted parallel to a 
given outgoing hard particle, in the case that this is massless. We will not be concerned with these 
here, but let us briefly remark that Ji is called a jet function, and contains collinear singularities 
associated with outgoing particle i. Because gluons can be soft and collinear at the same time, 
one double counts these singularities by including them both in the soft function S and the jet 
functions {Ji}- One corrects for this double counting by dividing by eikonal jet functions Ji, which 
can be thought of as the jet functions Ji evaluated in the soft limit. 


From now on we focus exclusively on the soft function S which, as described above, collects all 
singularities arising from multiple gluon soft gluon emission, including overlapping soft / collinear 
divergences. We will see that the soft function has a Feynman diagram intepretation by itself, 
so that it will be possible to calculate the soft function independently of the particular details 
of a given hard scattering process (other than the number of hard outgoing particles, and their 
momenta). Furthermore, we will see that the soft function S has an exponential form, such that 
one may schematically write 


S ~ exp 


Y^a^Cn 


( 3 ) 


n 


where a is the coupling constant, and the {cn} are coefficients which depend upon the momenta of 
the hard outgoing particles. This formula is a highly useful result, for the following reason: because 
soft gluon corrections appear in an exponent, they are summed up to all orders in perturbation 
theory. Expanding the exponential as a Taylor series in a gives terms to infinite order, where 
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IR divergent contributions at higher orders are related to singularities at lower orders. Successive 
terms in a in the exponent sum up successive towers of divergent contributions in the perturbation 
series. We see, then, that exponentiation is crucial to the resummation of soft gluon contributions, 
and that the problem of classifying IR divergences to all orders in perturbation theory is equivalent 
to the study of the exponent of eq. Q. 

We can go further than this. It turns out that, in both abelian and non-abelian gauge theories 
(as well as perturbative quantum gravity), the exponent of S can be given a Feynman diagram 
interpretation directly. That is, the cofhcients Cn can be obtained from a special class of Feynman 
diagrams. The particular properties these diagrams have depend on which theory we are consider¬ 
ing. For now, we will refer to any such diagram as a web: 


A web is a diagram that contributes to the exponent of the soft function. 


We will make this definition more precise (and even revise it) in the various theories to be examined 
in what follows, but this first definition of a web efficiently encapsulates the main idea. Before 
examining specific theories, it is useful to motivate why webs are useful in the first place. There 
are a number of reasons, which include the following: 


1. By definition, webs allow one to compute the exponents of soft gluon amplitudes directly. 
This is in contrast to the traditional way of calculating the exponent at a given order, which 
is to subtract from the (unexponentiated) amplitude those terms which arise from the expo¬ 
nentiation of lower order terms. This involves a large amount of unnecessary work, and it is 
much more efficient to be able to work directly in the exponent. 

2. Webs are, in an important sense, complete Feynman diagrams. Thus, some finite contributions 
may exponentiate as well as parts which are formally singular. Such finite contributions can 
be numerically large in many scattering processes, so that resummation of them is a useful 
means of increasing the precision of theoretical predictions. 


3. Recent work has begun to generalise factorisation formulae beyond the pure soft gluon 
approximation (also referred to as the eikonal approximation for the outgoing hard parti¬ 
cles) , 39 - ^ . Webs may be an essential ingredient of this programme in the most general 


case. 


4. Webs, and their associated mathematical structures, may give new insights into the structure 
of IR singularities in a variety of theories, including in particular QCD and J\f = A Super- 
Yang-Mills theory. 


5. There are links between webs in non-abelian gauge theories and interesting problems in the 
field of enumerative combinatorics. Some of these problems may have practical applications in 
computer science, thus there is potentially an intriguing possibility that lessons from quarks 
and gluons can be applied to e.g. network transfer protocols, or search algorithms. 


6. Recent progress hints at systematic methods for computing web diagrams 55 -64 . This offers 


the hope of being able to compute all-order properties in the exponents of amplitudes, which 
constitutes an unprecented look at intricate structures in perturbation theory. 
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Figure 2: Amplitude for an offshell photon decaying to a fermion-antifermion pair: (a) leading 
order; (b) dressed with a single photon emission. 


Having introduced the basic idea of webs, and also having hopefully conveyed why they might be 
useful, let us now begin to examine the simplest case of webs - namely, the exponentiation of soft 
photon corrections in abelian gauge theory. 

2 Abelian Exponentiation 

In this section we consider QED, and show in a specific example that soft photon contributions ex¬ 
ponentiate, where the exponent can be interpreted directly in terms of a type of Feynman diagram. 
This analysis is similar to that originally carried out in [65] . 

We begin by considering the simple scattering amplitude shown in[^a), consisting of an (off-shell) 
photon decaying into a massless fermion / anti-fermion pair. We label the 4-momenta of the latter 
by p and p respectively, as shown in the figure. Ignoring coupling factors for brevity, the amplitude 
in momentum space is given by 

A = u{p)-i>^v{p), (4) 

where p is the Lorentz index of the offshell photon. Let us now consider adding a single photon 
emission, as shown in figure [^b). The corresponding amplitude is 

(5) 

where Dai 3 {k) denotes the photon propagator, and we do not show explicitly the integral over the 
gluon momentum k. One may now consider the soft limit —?■ 0. In this limit, we can ignore the 
factors of in the numerator of eq. ([^, and neglect the terms in the denominators. Using also 
the fact that = 0, one finds 

Ai ( 6 ) 

We can simplify this expression by using the anticommutation property of the Dirac matrices, 
{7^)7^} = to write 

m(p) 7 " ^ = —u{p) jiy" -|- 2p°‘u{p) 
p^^v{p) = — 7 ^ pv{p) + 2p^v{p). 
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Figure 3: Pictorial depiction of the factorisation of the one loop graph of figure [^b) in the soft 
photon limit, where • denotes an eikonal Feynman rule. 


The Dirac equations for the fermion and antifermion spinors {u{p) ^ =pv{p) = 0) imply that the 
first terms of eqs. 0 and Q vanish. Then eq. Q becomes 

-Ai ^ [u{p)Yv{p)] (9) 

We may note the following from this expression: 

• The tree-level interaction factors out, and is not IR divergent. Thus, this is an explicit 
example of the factorisation of soft physics from hard, which is embodied more generally 
by eq. (§. The essential reason for this factorisation is that a soft photon has an infinite 
Compton wavelength, and therefore cannot resolve the details of the hard interaction (which, 
by the uncertainty principle, takes place over relatively short distance scales). 


The emission of a soft photon from the fermion or anti-fermion line is described by the eikonal 
Feynman rule 


pU 

p ■ k' 


( 10 ) 


where the sign is different on the two lines due to the fact that the momentum k is outgoing 
at one vertex, but incoming at the other. Here we have considered soft photon emission from 
fermion lines, but in fact the same Feynman rule would have resulted had we considered 
scalar or vector particles. The reason for this is again due to the large Compton wavelength 
of the soft photon, which cannot resolve the magnetic moment of the emitting particle, hence 
is insensitive to the spin. 


• We can think of the factors multiplying the tree-level interaction as a Feynman diagram 
by itself, albeit with effective Feynman rules. This is illustrated pictorially in figure]^ where 
we draw the one loop amplitude in the soft limit as the product of the tree level diagram 
with an effective Feynman diagram involving two eikonal Feynman rules joined by a photon 
propagator. The latter part is a subdiagram of the full amplitude which spans the hard 
external lines, and generates the IR singularities of the amplitude following the integration 
over the photon momentum k. 


Having observed the factorisation of soft photon physics at the one loop level, the next step is to 
consider what happens when we add any number of soft photon emissions. At a given order 
there are many different possibilities, two of which are shown in figure Qa). Each diagram gives 
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(a) 


(b) 



Figure 4: Two diagrams which contribute at where in the left-hand diagram we label the 

momentum of the gluon by ki. 



Figure 5: Pictorial representation of the sum of all soft photon diagrams at 0{a^), using the same 
notation as figure 


a somewhat complicated expression, as we will see in the following. However, a very simple result 
is obtained in the soft limit for the sum of all diagrams at 0(a"'), which is illustrated pictorially 
in figure This shows that in the sum of all soft photon diagrams, the tree-level interaction 
again factors out. The remainder is simply related to the one loop diagram raised to the power n, 
where the product is meant in the same sense as in figure]^ We may prove this result as follows. 

Let us first consider the particular diagram shown in figure [^a), consisting of n uncrossed photon 
exchanges. This is usually referred to as a ladder diagram in the literature. The expression for this 
diagram has the general form 

4"^ = (fci)... {kn), (11) 

where the tensors and collect all factors from the fermion and anti-fermion lines 

respectively, and we have included the product of photon propagators such that the photon has 
Lorentz indices a*, (3i. Again, we do not explicitly show the integrals over the loop momenta ki, 
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and ignore coupling constants etc. The factors associated with the quark line are 


which in the soft limit of ki — 




{p - kiY{p -ki- k2Y ...{p-ki-k2-...- knY 

• 0 for each i becomes 

u{p)'y‘^^ |i 7“2 _ . 7 “" |i 

[-2p ■ ki] [-2p ■ {ki + k2)] [-2p ■ {h + k2 + ... + kn)]' 


( 12 ) 


(13) 


As in the case of one photon emission, one may anticommute the factors of ^ through the Dirac 
matrices, to act on the spinor u{p). All terms involving u{p) p vanish from the Dirac equation, 
leaving 


p • A:i p • (fei + ^2) p • (fei + ... + A:„) ■ 
Similarly for the factors on the anti-fermion line one finds 


Q/^i'"^"n(p) = — 


pPi 


. p^ 


p ■ kip ■ {ki + ^ 2 ) ■ ■ - P ■ {k^ + 
such that the amplitude of eq. 0 becomes 

4 "> = a{- 


+ kn) 


v(p), 


(14) 


(15) 


^ p^n pl^\ 


pi^i^ 


~p-kip- {ki + k2). ..p- {ki + 
X Da\l3\{ki) . . . D„^,g^(A:„)|’. 


+ kn)p ■ kip ■ {ki + k2).. .p ■ {ki + ■ ■ ■ + kn) 


(16) 


Again we see that the tree-level interaction Aq has factored ont, and clearly this will happen for 
any such soft photon diagram at arbitrary order. Note that the denominators in eq. (16) indicate a 
complicated dependence on the photon momenta ki. These are coupled nontrivially along both the 
fermion and the anti-fermion line, so a given photon somehow knows about all other photons lying 
between itself and the final state hard particles. However, we have only considered one diagram at 
this order, and must now consider summing over all such diagrams. 

This set of diagrams is illustrated for the n = 3 case in figure where we see that there are six 
diagrams involving three photon emissions. We can think of obtaining this set by fixing the order 
of the photon emissions on the lower (anti-fermion) line, and summing over all permutations of 
emission orders on the upper (fermion) line. Thns, the total number of diagrams should be the 
same as the total number of permutations of three objects, which is 3! = 6 as observed. Returning 
to the general case of n photon emissions, we may write the complete sum over all diagrams from 
eq. (16) as 


.47 = 


tp • Kl p • 

xE 


(-!)> 




p^n p^^ 


.P^ 


{ki + k2)...p - {ki + ■■■ + kn) 

1 


p ‘ p * “1“ kT^2) • ‘ • P ‘ {k^Ti T ... T k^Y ^) 


k^aiPi (kl) ■ ■ ■ {kn) 

}■ 


( 17 ) 













Figure 6: Complete set of diagrams involving three photon emissions at 0{a^). 


Here vr labels a permutation of (1, 2,..., n), such that this maps to (vri, 7r2,..., Tr^). The rather 
formidable looking formula of eq. can be drastically simplified if we make use of the eikonal 
identity 


E 


1 

P * ^TTl P * (^TTl T ^TT2 ) ■ ■ • P * (^TTl T ■ ■ ■ “t“ ^TTri ) 



(18) 


In order to make sense of this result, it is instructive to consider the case of n 
reduces to 


2, in which eq. (18) 


1 

p- kip- {ki + k2) 


+ 


1 

p- k2p- {ki + k2) 


p ■ ki + p ■ k2 
p ■ kip ■ k2P ■ {ki + k2) 


1 

p ■ kip ■ k2' 


(19) 


which is indeed the stated result. We see that the eikonal identity is a sort of partial fraction 
identity, and the general result could be proven e.g. by induction from lower orders. 


Problem 1. Show by explicit calculation that the eikonal identity works 
for the n = 3 case. 


Substituting eq. (18) into eq. ([Xt]), the latter becomes 




(- 1 )>° 


.p 




.P^ 


p ■ kip ■ {ki + k2). ■ ■ p ■ {ki + ... + kn) 


Daifii (kl) ■ . . Donfin (^n 



(20) 
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This result is already a lot simpler than eq. for one important reason: the photons on the upper 
line are now completely decoupled from each other. This is not true for the photons on the lower 
line. However, one may achieve this by exploiting the fact that all of the momenta ki are dummy 
variables which are ultimately integrated over. Furthermore, the product of photon propagators 
contracted with the momentum factors is symmetric under interchange of any two momenta, subject 
to suitable relabelling of dummy Lorentz indices. Thus, by relabelling the integration variables one 


may replace the particular combination of photon momenta appearing in the first line of eq. (20) 
with any permutation of these momenta. That is, one may write 


(fki. 




1 




p ■ kip ■ {ki + k2).. .p ■ {ki + ... + kn) 

1 


P ' kjYi p ‘ {k-Ki T ^772) ‘ • 'P ' (^TTi T .. . T 


( 21 ) 


In the second line, we have used the fact that each term gives the same result (by relabelling 
dummy integration variables), and that there are nl such permutations. Substituting this result 


into eq. (20) and tidying things up (again omitting the integrations over the momenta for brevity), 
one ultimately finds 




A-n 


\P-hj \P- k. 


P^i 


(ki) 




nl 


-p 


p■k J \p■ k 


pl3 


Dap{k) 


( 22 ) 


which is precisely the result expressed pictorially in figure namely that the sum over all n photon 
graphs in the soft limit is simply related to the power of the one loop graph. 


The astute reader will have noticed that eq. ( p2) ) (which includes the leading order result) is the 
term in the Taylor expansion of an exponential. Thus, one may explicitly compute the sum 
of all diagrams involving any number of single soft photon emissions, to all orders in perturbation 
theory! The result is 


^ = ^0 exp 


dfik f —p°‘ 

{2'kY VP • k 



p ■ k 




(23) 


where we have reinstated the integral over the gluon momentum. This result is shown using our 
pictorial notation in figure In later sections, we will end up considering diagrams with many 
hard particles connected by soft gauge bosons. Thus, it is useful to redraw figure so as to include 
the hard interaction in the exponent, as shown in figure The purpose of this is simply to clearly 
display which lines the photon is emitted between - we do not mean that the hard interaction itself 
exponentiates. 


We have now derived our first example of exponentiation. Some further comments: 


(i) The one loop graph in the exponent of eq. (23) has an IR divergence. Were we to explicitly 
expand the exponential, we would thus find that IR singularities at higher orders in the 
amplitude are related to singularities at lower orders. 
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Figure 7: Pictorial representation of the exponentiation result of eq. (23). 



Figure 8: Alternative drawing of figure [7j making clear which outgoing particles emit the soft 
photon. 


(ii) Here we considered only graphs with single photon emissions, finding that these are entirely 
given by the exponential of the one loop graph. In fact, this is not the whole story - soft 
photons could also connect off the external lines, via fermion bubbles. The complete result is 
then that the exponent contains all connected soft photon subgraphs, as depicted in figure]^ 
In principle one can prove this using a generalisation of the above method [^. We will 
instead prove this result using alternative means in the following section. 

(iii) Above, we only considered two outgoing hard particles. However, the result can be extended 
to any number of external lines, where again the exponent contains all possible connected 
subgraphs which span the hard particle lines. 

(iv) Things are, unsurprisingly, more complicated in non-abelian gauge theories. Crucial to the 
above argument was that we could sum over all permutations of photon emissions, where these 
were weighted equally. This is no longer the case in e.g. QCD, where the emission vertices 
carry colour matrices T“ which do not commute. Thus, an additional level of combinatoric 
complexity arises. 


(v) 


Things are not much more complicated in perturbative quantum gravity. This is a theory 
that shares features of QED and QCD. It is QED-like in that soft graviton emission vertices 
commute with each other (i.e. do not carry a non-commuting matrix-valued charge). It is 
QCD-like in that there are three and four graviton vertices, which mimic the correspond¬ 
ing gluon self-interactions in QCD. However, the charge in gravity is the 4-momentum of 
a particle. Thus, all multiple vertices between soft (zero momentum) gravitons vanish. A 
careful analysis also shows that diagrams in which gravitons couple to matter loops also van¬ 
ish p7 29 . Hence, IR divergences are given to all orders in quantum gravity purely from 


the exponentiation of the one-loop result, a property known as one-loop exactness. In other 
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Figure 9: Complete result for abelian exponentiation, in which all possible connected soft photon 
subgraphs appear in the exponent. 


words, the structure of IR divergences in quantum gravity is, remarkably, simpler than in 
QED! 

In section we defined a web as being a diagram which contributes to the exponent of the soft 
function. We have now seen our first examples of webs, and the results of this section can be 
compactly summarised as follows: 


Webs in abelian gauge theories are connected subdiagrams. 


At this point, we could leave the case of QED in favour of exploring non-abelian gauge theories. 
The traditional treatment of webs in the latter case proceeds by generalising the eikonal identity 
to take into account non-trivial colour structure. One may then explicitly solve for the structure 
of webs in scattering amplitudes which involve only two coloured external particles 66-68 ^ We 
choose not to follow such an approach here. Firstly, there is the lack of generality implied by the 
restriction to processes involving only two coloured particles (although three particles can also be 
treated straightforwardly using similar methods Secondly, one may develop a more general 

approach for deriving exponentiation, which works for both abelian and non-abelian gauge theories, 
and any number of outgoing particles. In order to introduce this method, we will start by rederiving 
the results of abelian exponentiation. This is the subject of the following section. 


3 The Path Integral Approach 


In the previous section, we saw that the soft function itself has a Feynman diagram interpretation, 
where these diagrams consisted of effective Feynman rules for emission of soft gluons from hard 
outgoing particles. The propagators joining these vertices were the usual propagators arising from 
the action for the soft gauge field. Furthermore, these diagrams were subdiagrams in the full am¬ 
plitude. That is, if one removes the hard interaction and the hard outgoing particle lines, what is 
left is the soft gluon diagram. This is illustrated in figure [Tol for a scattering amplitude involving 
five external particle legs. In practice, as we saw at the end of the previous section, it is useful to 
also draw the external lines in the soft photon subdiagram, so as to indicate which external line a 
given soft photon was emitted from. This matters given that the hard momentum enters the 
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For a pe daRORic al exposition, see 


See also 


70 


72 


69 


for an examination of the general multileg case. 
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(b) 


(a) 

Figure 10: (a) An example scattering amplitude involving five external hard particles, and two soft 
photons; (b) soft photon subdiagram. 


eikonal Feynman rule. 


Given that Feynman diagrams exist purely for the soft subdiagrams, it follows that there must 
be a quantum field theory for the soft gauge field by itself, whose diagrams are the soft photon 
subdiagrams. One way of specifying a QFT is to write down its generating functional, from which 
all Feynman diagrams can be derived. Here we will simply write down the generating functional 
which generates soft photon subdiagrams (see 12 for a full derivation). It is given by[^ 


Z = / VA,, 


JS[A^ 


H' 

k=l 


J J dx'^Af,(x,,) 


(24) 


where is the (soft) gauge field with action and the product is over external particles, 

L in total. The exponential factors contain the space-time trajectory of the particle, which is 
given by 

< = Kt, (25) 

where is a (constant) 4-vector proportional to the 4-momentum of the particle, p^. Note that 
we have taken each outgoing particle as having a zero initial position. Thus, all the hard particles 
originate from the origin, which corresponds to having shrunk the hard interaction which produces 
these particles to an infinitely small region of space-time. This is an approximation of course, but 
it turns out that corrections to this approximation are formally subleading, in that they involve 
additional powers of the emitted gluon momenta, which in the soft limit have been taken to zero. 
The physics of this approximation is that soft gluons, as discussed above, have an infinite Compton 
wavelength. Thus, they are not able to resolve the internal details of the hard interaction, which 
consequently appears pointlike. 


Let us discuss the form of the generating functional of eq. (24). Firstly, there is a path integral 
over all possible configurations of the soft gauge field, where each is weighted by a phase factor 
involving the classical action. This is what one expects for the path integral representation of a 
quantum field theory. The additional factors (one for each external leg) are exponential, where the 
exponent is linear in the gauge field. This thus acts as a source term for the soft gauge field. In 
other words, vertices will be generated, which enable a soft photon to pop out of the vacuum. The 


^Here we have absorbed a factor of the coupling constant into the gauge field. 
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form of the exponent tells us that each vertex is located in space-time along one of the external leg 
trajectories x^. We can see what the momentum space Feynman rule is for such a vertex by first 
rewriting the exponent of each external line factor as 


exp 


dx^A^ 


= exp 


poo 

i / din'll A fj,{x^) 

. Jo 


(26) 


where we have substituted the trajectory of eq. (25). Next, we may write the soft gauge field in 
terms of its Fourier transform using 


= 


d^k 

W) 


^A^{k)e 


ik-x 


(27) 


Substituting this into eq. (28), the exponent becomes 

d^k 


i / dtn1,A^{n'^t) = 


{2'kY 

d'^k 
d<^k 


AYk)i / 

-| oo 


A^{k)i'n}t 


ik ■ Uk 


J 0 


Mk) - 


nl 


Uk ■ k 


(28) 


J {2'kY ^ 

Here we have been rather cavalier in dealing with the upper limit in the second line. A more 
careful treatment - including the Feynman ie prescription - gives the same result. By the usual 
procedure for reading Feynman rules from the exponent inside the path integral, the integrand of 
the momentum integral (represented by the term in brackets in the last line of eq. (28)) is the 
momentum space Feynman rule. Given that this is evidently invariant under rescalings of the 
vector oc p^, we may replace this Feynman rule with|^ 

Pk 


Pk-k' 

which is the correct Feynman rule for the emission of a soft photon from a hard external particle, 
as given in eq. (10). 


Problem 2. For hard particles emitting soft gravitons (described by a spin 
two field /i^[y(x)), the analogue of the external line factor of eq. (26) is|^ 


exp 


K I 

-i- dsp^^p'" h^Ysp) 

^ Jo 


where p^ is the momentum of the particle emitting the gravitons, s is a 
parameter such that x^ = and k = -v/IhvrGjv, where Gat is Newton’s 
constant. Show that this leads to an effective momentum-space Feynman 
rule for graviton emission given by 

np^ p’^ 

2 p ■ k 


5 

16 


The rescaling symmetry of the eikonal Feynman rules is actually broken when the Wilson lines are null - see e.g. 
for a discussion. For our purposes, we can take the Wilson lines to be massive i.e. p\ yi; 0. 
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The above discussion justifies why eq. (24) is the correct generating functional for soft photon 
snbdiagrams. It generates the right Feynman rnles for emission of soft photons from external lines. 
Fnrthermore, these will be joined by propagators, as generated by the action for the soft gange 
field Thus, the diagrams generated by eq. (24) are indeed the soft photon diagrams we 


are after, which are in turn snbdiagrams in the fnll amplitude. 


There is another way to look at the external line factors in eq. (24), which is to notice that they 
are Wilson line operators. The general definition of such an operator is|^ 


<I>(C) = exp 


i / dx^ A^{x) 


(29) 


where C is some contour in space-time. Such operators transport gauge information from point to 
point in space-time, along the chosen path. What eq. (24) then generates is a vacunm expecta¬ 


tion valne of Wilson line operators, evalnated along the classical trajectories of the outgoing hard 
particles. This is in fact a restatement of an old result [5,75 , namely that IR singularities due to 


soft gluon emission from hard qnarks or gluons can be described by dressing the hard particles by 
Wilson line operators. The physics behind this is as follows: if outgoing particles are emitting soft 
(zero-momentum) gluons, they do not recoil. Thus, they must follow their classical straight-line 
trajectories. They can, however, change by a phase. If this phase is to have the right gauge- 
transformation properties to be part of a scattering amplitude, it must be a Wilson line operator. 


Given that eq. (24) generates the soft function in the amplitude, another way to say the above is 


that the soft function is a vacuum expectation value of Wilson lines, or: 


II'^‘ 

k=l 


(30) 


where is a Wilson line operator evaluated along the trajectory Xk of particle k. Indeed, this is 
usually taken as the definition of the soft function (see e.g. (I^). 


Having introduced the generating functional for the soft gauge field theory, we now proceed to 
use this in order to rederive abelian exponentiation. To this aim, we introduce a crafty statistical 
physics approach, which is the subject of the following section. 


3.1 The replica trick 

In the previous section, we introduced a generating functional for a QFT for the soft gauge field 
A^, which generated soft photon subdiagrams. If we consider the simple case of QED with no prop¬ 
agating fermions (so that soft photon graphs contain no fermion bubbles), the path integral over 


An in eq. (24) is Gaussian, and can be carried out exactly. This would lead immediately to abelian 


exponentiation, which we already enconntered in section In the more complicated cases of QCD 
and QED with fermion bubbles, it is no longer the case that the path integral in the generating 
functional can be performed exactly, and so we instead wish to introduce a method for deriving 
exponentiation of soft gauge boson diagrams which will be fully general. Specifically, we will use 


^Note that we have been somewhat lapse with our notation for this action, which also includes couplings to matter 
particles, so that fermion loops are generated in soft photon diagrams. 


8 


As in the preceding discussion, we have here omitted a factor of the charge e of the emitting particle. 
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the replica trick, a technique borrowed from statistical physics (see e.g. [^). We will see that this 
provides an elegant shortcut to exponentiation, which bypasses the nontrivial combinatorics (e.g. 
the eikonal identity) that we encountered in section]^ 


The argument proceeds as follows. Instead of the field theory defined by eq. (24), let us instead 
consider a theory with N identical copies, or replicas, of the soft gauge field A^. We can label these 


replicas by \ where 1 < j < N. Importantly, we take these replicas to be non-interacting. That 
is, a given replica does not interact with other replicas, so that the set of replica gauge fields is 
mutually independent. 


By generalising eq. (24), the generating functional for the replicated theory is 


Zrep, = exp 


N 




f=i 


Lfc=l 


J f dxk-AW 


Lfc=l 


J f dxk-AW 


(31) 


Here there is a path integral over all possible field configurations of each replica. The total action 
for the replicas is the sum of all the individual replica actions, as follows from the fact that the 
replicas are non-interacting. Furthermore, we now have a product of Wilson line factors for each 
external line, where this product occurs for each replica number from 1 to N. We can simplify 


eq. (31) by noting that the Wilson line factors in an abelian theory commute with each other, so 


that they can be combined into a single exponential to give 


Z,ep, = VAWvA(^K..VAI^^ exp 


N 




i=i 


n 

k=l 


exp 


N 


i ^ / dxk ■ A^^'l 




(32) 


What do the diagrams generated by this theory look like? Ignoring the replica index, the generat¬ 
ing functional looks similar to the non-replicated theory, and so diagrams in the replicated theory 
have similar topologies to the original soft photon theory. However, we have to take account of the 
fact that each connected piece of a given diagram could have a different replica number. This is 
illustrated in figure [TH which shows three examples of soft photon diagrams, dressing the hard in¬ 
teraction of figure j^a). The first example consists of two single soft photon emissions, each of which 
may potentially correspond to a different replica (shown using different colours in the figure). We 
can also think about this in terms of assigning indices i and j to each photon, where 1 < i, j < N. 
In the second example, there are two disconnected pieces in the soft photon subdiagram: a single 
photon emission, and four photons connected to a fermion bubble. Each of the separate connected 
pieces is constrained to have a single replica number. This follows from the fact that the action for 
each replica, implicitly includes coupling to fermion bubbles etc. Thus, replicas cannot mix 

via matter bubbles, which would ruin the fact that the replicas are independent. The third example 
in the figure shows two similar connected pieces arranged in a different way, so that there is some 
sort of non-trivial kinematic entanglement between the different parts of the photon subgraph. This 
does not affect the fact, however, that the two connected pieces are completely separate from each 
other if the external lines are removed (as defines the soft photon subgraph), and thus they may 
still have potentially different replica numbers. 


The above discussion is perhaps rather unmotivated and mysterious! Let us now see how we can use 
the replicated theory to prove abelian exponentiation. Firstly, note that the replicated generating 
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(b) 


(c) 




(3) 


Figure 11: Examples of soft photon diagrams in the replicated theory of eq. (31). Different colours 
represent potentially different replica numbers. 


functional of eq. (31) is related to the generating functional of the original theory, eq. (24), by the 
relation 


Zrep. — Z 


N 


(33) 


where as usual N is the total number of replicas. This follows from the fact that the replicas are 
non-interacting, so that one could choose to rewrite eq. (24) as a product of path integrals over 
each replica individually, thus making the structure of eq. (33) manifest. We may now Taylor 
expand this equation in N to give 


Zrep = l + NlogZ + 0{N^). 


(34) 


The left-hand side of this equation corresponds to a complete sum of Feynman diagrams in the 
replicated theory. These will, in general, be N dependent. Equation (34) then tells us that if we 


expand each Feynman diagram as a function of N, the coefficient of the 0{N) piece contributes 
to the logarithm of Z (the generating functional of the original theory). Note that this argument 
does not depend on N being small - higher order terms in N contribute to higher powers of the 
logarithm of Z, which is irrelevant to our present logic. 


It follows from the above discussion that if we find the 0{N) part of all diagrams in the replicated 
theory (let us call this W for a given diagram), then the original generating functional is given by 


Z = exp 


. w 


(35) 


Thus, Z has a manifestly exponential form. This is not a particularly impressive statement in 
general - it is always true that Z is the exponential of its own logarithm! However, what the 
replica trick gives us is an explicit procedure for determining the logarithm of Z; in other words, 
for determining which diagrams exponentiate in the non-replicated theory. 


To see how this works in this abelian case, consider for example the diagrams in figure I13b) 
and[IT|(c). These both have two separate connected pieces in the soft photon part of the diagram, 
and each of these has a potentially different replica number. Given that there are N choices for the 
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first replica number and a further N for the second, it follows that diagrams with these topologies 
are 0{N‘^). Generalising this argument, a diagram with m separated connected pieces is 0{N'^). 
Such diagrams, by the above argument, do not contribute to the exponent of Z, as they have no 
part which is linear in N. The only diagrams which are linear in N have m = 1 i.e. have only one 
connected piece. We have thus derived, from the replica trick, that connected diagrams exponen¬ 
tiate in the soft photon theory All disconnected diagrams, such as those shown in figure [TT|^b) 
and (c), are obtained by the explicit exponentiation of lower order diagrams. 

Recalling that Z generates subdiagrams in the full amplitude, namely those parts of the graph that 
contribute to the soft function, we have derived the result 


S =exp [Connected soft photon subgraphs] . 

This is the abelian exponentiation result that we (partially) derived using the eikonal identity in 
the previous section. Some further comments are as follows: 


(i) The replica trick proof is much simpler than the eikonal identity approach. Essentially, what 
the replica trick does is to hide the tricky combinatorics associated with building up the ex¬ 
ponential series. Using a traditional combinatoric approach would be even more complicated 
in non-Abelian theories, although it can of course be done 1^-68 70 77|. 


(ii) The replica trick works for any number of external lines. This will be particuarly useful in 
the case of non-abelian gauge theories, in which results have traditionally been separated 
into the case of amplitudes with two coloured particles, and the general case of multiparton 
scattering. 


(iii) The above result shows that all connected diagrams exponentiate, including those containing 
fermion bubbles. We would have had to have proved this separately in the eikonal identity 
approach of section (although admittedly this is a minor modification [65|). 

(iv) In the replica trick analysis, we explicitly used the fact that the Wilson line operators commute 
with each other, in order to combine their exponents into a single exponential factor for each 
external line. A similar result holds in quantum gravity, but not in non-abelian theories such 
as QCD. The latter have colour matrices associated with the emission of a soft gluon, which 
complicates the combinatorics of the path integral. In particular, the standard textbook 
result of exponentiation of connected diagrams does not hold for the field theory describing 
a soft gluon field (as we will see). 


Having introduced the replica trick in this section, and described how this can be used to derive 
exponentiation properties, we are now ready to explore exponentiation in non-abelian theories. 


4 Non-abelian exponentiation 

In the previous section, we introduced the replica trick and showed how this can be used to de¬ 
rive abelian exponentiation. The main idea was that the 0{N) part of Feynman diagrams in the 

®This is in fact a special case of a general result in quantum field theory (which can be found in many textbooks, 
albeit often in an appendix), that connected diagrams exponentiate. 
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Figure 12: (a) Schematic depiction of qq —>■ qq scattering, where the labels denote colour indices; 
(b) The two independent ways of joining up the colour flow. 


replicated theory contributes to the exponent of the soft function in the non-replicated theory. 
Thus, one has an explicit procedure for calculating the exponent of the soft function directly. In 
non-abelian theories, the same overall procedure can be used, but leads to a much richer structure 
in the exponent that is still not fully understood. 

Before examining non-abelian exponentiation in detail, it is worth making a few remarks regarding 
the existing literature on this subject. These are given in the following section, which can be skipped 
if necessary without affecting the flow of the argument from the previous section. However, what 
follows is useful for reasons of completeness, and for making contact with previous work on webs. 


4.1 Two parton vs. multiparton scattering 


Historically, non-abelian exponentiation was first considered in the case of scattering amplitudes 


involving only two coloured particles 66-^ (plus any number of colour singlet particles). Physical 


examples of these processes include Drell-Yan production of a vector boson, Higgs boson production 
via gluon fusion, deep inelastic scattering (in which an electron collides with a proton, which 
subsequently breaks up), and electron-positron annihilation into a quark-antiquark pair. In such 
processes, each external quark or gluon has a colour index. The special nature of amplitudes with 
only two coloured particles is that there is only one possible colour flow at the hard interaction 
vertex: the colour indices of the two quarks or gluons involved must join up with each other (by 
colour conservation). In a general multiparton scattering process, there are many colour indices 
joining up at the hard interaction vertex, and consequently there are many different ways of joining 
them up. Consider, for example, the case of qq —)• qq scattering, as shown in figure [T^ a). Whatever 
is happening in the hard interaction, it remains true that there are only two ways of joining up 
the colour indices of the incoming and outgoing quarks. These are shown in figure [l2|(b). The 
scattering amplitude of this process will have four colour indices, as labelled in the figure [T^a). 
The two ways of joining the colour flow means that this amplitude must have the following form: 


■^ijkl — -^1 hk^jl + -42 hl^jk-, 


(36) 
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Figure 13: Examples of: (a) a two-eikonal line irreducible diagram; (b) a reducible diagram (the 
red line separates the two connected pieces). 


where 6ij is the Kronecker delta symbol, and {Ai} are kinematic coefficients which are independent 
of the colour index structure. We may rewrite this equation as 

A = '^Aici, (37) 

I 

where 7=1,2 and the colour tensors {c/} are defined via 


Cl — Sikdji, C2 — SiAjk- 


(38) 


In other scattering amplitudes, there are more ways of joining up the colour flow in general. How¬ 
ever, amplitudes still have the generic form of eq. (37), with an appropriate set of colour tensors. 
These are not uniquely defined - any basis obtained via a linear transformation will do. Having 
chosen a basis, we can work with the components of the amplitude Ai, and we see that amplitudes 
are vectors in colour-flow space. The soft function contains gluons which transfer colour from one 
external line to another in general. Thus, it mixes the various colour flows present in the hard 
interaction, and is matrix-valued in colour space. How such matrices act on the colour structure of 


the hard interaction is an interesting problem by itself - see e.g. refs. 78,79 


Having introduced the above terminology, we can see why soft gluon emission in the case of two 
hard external lines is simpler in general than the multiline case. In the two line case, there is only 
one colour tensor {6ij, where i and j are the colour indices of the two external lines). The soft 
function is then a 1 x 1 matrix, and different contributions to the soft function commute with each 
other. This is in contrast to the general case of multiparton scattering, in which the soft function is 
a non-commuting object in colour flow space, and has implications for the structure of the exponent 
of the soft function. It turns out that in the two-line case, one can classify which diagrams appear 
to all orders in perturbation theory using a simple topological criterion, namely that they are two- 
eikonal line irreducible. This means that one cannot disconnect a given diagram by drawing a single 
line through both external lines, and an example is given in figure [T^ Also shown is an example of 
a reducible diagram. In this case, figure [T^b) would not enter the exponent of the soft function, 
whereas [T^ a) would. Recall that diagrams such as the latter were not present in the exponent in the 
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abelian case, so that even in the two-line case QCD is already much more complicated than QEDp^ 


Problem 3. Draw all irreducible diagrams at 3 loop order. 


The above results for the two line case were derived in 166 -68 , and indeed it was a particular 


way of drawing two-eikonal line irreducible diagrams in those papers that first led to the term 
web. We will not present the derivation of this result here. It proceeds via a generalisation of the 


eikonal identity of eq. (18), which takes account of the fact that successive gluon emissions are 
non-commuting (for a recent presentation, see [^). Rather, our aim is to present a formalism for 
non-abelian exponentiation, which is set up for the general case of multiparton scattering. The 
special properties of two-line scattering will then emerge as a special case. The following sections 
are based heavily on 77 , but note that multiline webs have also been considered from a different 
point of view in 


For a discussion of how these two approaches are related, see appendix A 
of [80|. More recently, another approach has also appeared |71[|72|. 


4.2 The replica trick in a non-abelian theory 

We now return to elucidating the structure of non-abelian exponentiation. In the abelian case of sec¬ 
tion we used the replica trick to elegantly derive the exponentiation structure of the soft function, 
namely that connected soft photon subdiagrams enter the exponent. This strongly suggests that we 
should use the replica trick also in the non-abelian case. Let us then try this and see what happens. 


First we need to write down the generating functional for the soft gluon theory, whose diagrams 
are the soft gluon subdiagrams of a non-abelian scattering amplitude. This is a straightforward 


generalisation of eq. (24), and the result is 


^ = / VA,. 




UVe 


if dx'^A'^ 


(39) 


k=l 


This has a path integral over the soft gluon field, with a weighting factor containing the classical 
action. For ease of notation we use the same symbol S[Afj] as in the QFD case. However, the 
non-abelian action will contain three and four gluon vertices etc. in addition to the propagator and 
couplings to matter bubbles. Each external line, as before, has a Wilson line operator associated 
with it. This is similar to the QFD Wilson line, except for the fact that the gauge field A^ = 
is matrix-valued, where T^ is a colour matrix acting on the partonic colour indices of the line. 


We do not write these colour indices explicitly in eq. (39), but we do add an index k to the soft 
gauge field to remind us that this lives in the colour index space of the k^^ parton line. Given 
that the exponent of the Wilson line operator is matrix-valued, the exponential itself is ambiguous. 
Upon carrying out a Taylor expansion, it is not clear how to order the colour matrices in the correct 


fashion. This is corrected by the operator V in eq. (39), which stands for path ordering. This opera¬ 


tor is defined such that, in a given term in the expansion of the exponential, the colour matrices are 
ordered along the line according to the sequential ordering of gluon emissions. Parametrising the 


contour of a given line according to eq. (25), for example, this means that the colour matrices will 


“This is only one complication of course. There are also soft gluon graphs in QCD which contain three and four 
gluon vertices, which do not exist in QED. 
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be ordered according to the order of the t values at which gluons are emitted. This is a somewhat 
subtle point that is absent in the case of an abelian gauge theory, where no non-commuting charges 
appear. 


Applying the replica trick, one may consider the theory 


Z,ep. = / PAW...PAWexp 


N 

i=i 


P exp (i I dx'^ 

.k=l 


llVexp(i I 
Lfc=l 


(40) 


where we have again taken N independent (mutually non-interacting) copies of the gauge field. 
At this point in the abelian case, we collected together all Wilson line exponents to make a single 
exponential. We cannot do that in this case, due to the fact that the exponents of the Wilson line 
factors do not commute, if these act on the same external line. We can, however, commute Wilson 
line operators which act on different external lines. These have completely different colour indices, 
so that matrices on different lines commute simply by not living in the same colour algebra. This 


allows us to gather together the contributions from each external line separately in eq. (40), to get 




rep. 


[ PAW...PAW exp 

Co 

' _ 1 

L 

n 

n^exp (i f dx'^A’^A 

J 


k=l 

j=i \ J 


(41) 


To recap, the first product in this formula is over all external lines. Then, each external line factor 
has a product of Wilson line operators acting on it, one for each replica. Furthermore, the replica 
number in this product increases from left to right, from 1 to N. 


We cannot immediately read off the Feynman diagrams for the soft gauge held from eq. (41), as this 


contains a product of path-ordered exponentials for each external line. Were we inclined to rewrite 
this as a single path-ordered exponential, we would then be able to analyse the diagrams in more de¬ 
tail, by analogy with the case of the original theory. In general combining a number of exponentials 
involving non-commuting objects involves multiple applications of the Baker-Campbell-Hausdorff 
formula (see e.g. refs. (7T||72p l for approaches to exponentiation that build on this idea). However, 
we can circumvent this in a crafty way by exploiting the fact, mentioned above, that the replica 
number increases sequentially from 1 to in the product of Wilson line operators associated with 
each external line. Thus, one may write this product as 


N 

Pexp 

i=i 


dx^ A^^^ 1 = nv 


exp 


N 


•E 

i=l • 


Af > 


(42) 


Here TZ \s a repliea-ordering operator. Given a product of colour matrices on a given line, it 

ensures that the replica number increases along the line. We can define its action more formally by 

(i) 

considering two colour matrices acting on line k. Denoting by T),^ a colour matrix acting on line 
k and associated with a gluon with replica number i, one has 




rp(d rp(i) 


'Ti(d lyd) 


i < j 

i > j 


(43) 
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That is, TZ does nothing if the replica numbers associated with colour matrices are already increas¬ 
ing (or the same) along a given line, but reorders them if this is not the case. It is straightforward 
to generalise this definition to higher numbers of matrices. The TZ operator then acts on each term 
in the expansion of the exponential on the right-hand side of eq. (42), and reorders the colour ma¬ 
trices such that the replica numbers in each term are in increasing order. Then the right-hand side 
agrees with the expansion of the multiple exponentials on the left-hand side which, if carried out, 
will have the replica numbers similarly ordered in every term. Note that the notation of eq. (42) 
implies that one first path orders (according to "P), and then replica orders (according to TZ). In 
cases where these orderings are different 7Z overrides V. 


Substituting eq. (|42|) into eq. (|41|), we may rewrite the latter as 

N 


^rep. — /p4i)...P4^)exp 






wn 

k=l 


Pexp 



(44) 


Here, to simplify notation, we have defined an overall replica-ordering operator 7Z which is under¬ 
stood to act appropriately on each line k. Note that this now looks more similar to the QED case 
of eq. ( |32[ ), apart from the colour-specific details (which are also present in the non-replicated QCD 
theory), and also the presence of the TZ operator. 


Consider now a particular soft gluon diagram D with a given topology. In the original (non- 
replicated) theory, this has the general form 

F[D)C{D), 

where C{D) is the colour factor of the graph (collecting all the soft gluon colour matrices on each 
external line), and F{D) contains all the kinematic dependence. In the replicated theory, the same 
diagram topology will have the similar form 

FiD)C{D), 

where F{D) is the same as in the original theory (as in the QED case, the kinematics of a given 
diagram does not care about the assignment of the replica numbers). However, the colour factor 
C{D) is not the same as in the original theory. Firstly, it depends on the number of replicas N. 
Secondly, colour factors are different in the replicated theory due to the fact that colour matrices, 
which obey path ordering only in the original theory, get reordered by the TZ operator. The details 
of this reordering depend on the particular assignment of replica numbers in a given graph. To find 
the total colour factor in the replicated theory, one must also sum over the total number of replica 
assignments. The whole procedure is perhaps best illustrated by example. 


Consider the two diagrams at two-loop order shown in figure 14 
amplitude with four external lines 


Here we have considered an 


and shrunk the hard interaction to a point. In the replicated 
theory, the two diagrams in the figure give contributions to the (replicated) amplitude 


Aa = F(a)C(a), A = F(b)C(b) 


(45) 


respectively. To find the total colour factor of diagram (a), one must consider all possible assign¬ 
ments of replica number i and j, as these will potentially contribute different amounts to the total. 
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Figure 14: Example two-loop diagrams, where i and j denote replica numbers. 


h j 

Replica-ordered colour factor 

Multiplicity 

i = j 

C{a) 

N 

i < j 

C{b) 

N{N -l)/2 

i > j 

C{a) 

nIn -l)/2 


Table 1: Replica trick analysis of the diagrams of figure [141 3.s described in the text. 


depending on the action of the replica ordering operator TZ. Given that the action of this operator 
only depends upon whether replica indices are greater or less than each other (or the same), we 
need only consider three cases: i = j, i < j and i > j- These are given in the first column of table 
and the third column gives the total number of replica number assignments corresponding to each 
hierarchy of i and j. In the second column, we note the A^-independent part of the contribution 
to the colour factor for each diagram. Considering first the case of z = j, this is simply the colour 
factor C{a) of (a) in the original theory, as no replica reordering takes place if the replica indices are 
the same. If z < j, the gluons are not ordered in terms of increasing replica number on the upper 
right-hand line (where, according to the arrows to designate our choice of momentum flow, we are 
taking the direction of the Wilson line as increasing away from the hard interaction vertex). Then 
the colour matrices get reordered, which produces the conventional colour factor C{b) of diagram 
(b). This is the reason why we have had to consider two diagrams at once - they mix with each 
other due to the TZ operator. Finally, we must consider the case i > j- In this case the replica 
ordering constraint is already satisfied, and thus the ^independent part of the contribution to the 
colour factor is C{a). 

Putting things together, the total colour factor in the replicated theory is 

C(a) = NCia) + [C{a) + C(6)] 

_ , r C(a)-C(6) 

[ 2 

Here we have formed the total by multiplying together the results of the second and third columns 
in table and summing these up. In the second line we have taken the 0{N) part of the colour 
factor, anticipating that this will be relevant in what follows. 


+ 0{N^ 


(46) 
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Figure 15: Diagrams considered in problem 4. 


One may carry out a similar analysis for diagram (b) in figure 14 and the total colour factor in 
the replicated theory is found to be 


C{b) = N 


C{b) - Cia) 


+ 0{N^ 


(47) 


We may now use the replica trick argument from section which said that the contribution of a 
given diagram D to the exponent of the generating functional in the non-replicated theory is given 
by the 0{N) part of this diagram in the replicated theory. Thus, we see that both diagrams (a) 
and (b) contribute to the exponent of the soft function in the original theory, and they do so in the 


combination (taking the coefficients of the 0{N) terms in eq. (46, 47)) 




Note that this can be rewritten in the matrix form 


Ha) 

m 


1 -1 

-1 1 


C{a) 

C{b) 


(48) 


(49) 


This is an interesting structure, as the vectors on the left and right-hand side (defined in the space 
of diagrams) involve properties of the original theory, as they should now that we have returned 
from the replicated theory. The matrix in the middle consists of constant numbers, and describes a 
mixing between the colour and kinematic degrees of freedom of soft gluon diagrams in the exponent 
of the soft function. This structure is in fact entirely general, as we describe in the following section. 

Problem 4. Without detailed calculation, write down the contribution to 
the soft gluon exponent of the diagrams in figure [T5| 


4.3 General structure of multiparton webs 

In the previous section, we examined a particular set of diagrams that mix under the action of 
the TZ operator. Both diagrams were found to contribute to the exponent of the soft function, but 
with a mixing matrix that entangles the colour and kinematic degrees of freedom. To see how this 
works in general, note first that the two diagrams of figure 14 are related by gluon permutations 
on the external lines. In this simple case, there is only one non-trivial permutation possible, on the 
upper right-hand line. Furthermore, this is a closed set: acting with more permutations will not 
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generate any additional diagrams. In general, given a diagram D involving arbitrary numbers of 
gluon attachments on arbitrary external lines (where these gluons may also be connected away from 
the external lines by fermion bubbles or multiple gluon vertices), we can always form such a closed 
set of diagrams related by gluon permutation operations. A further example is given in figure 
It is clear that a given diagram can only mix with diagrams from the closed set to which it belongs. 
This is because the source of the mixing is the TZ operator in the replicated theory, which is itself 
a form of gluon permutation operator, transforming between diagrams in such a closed set. 

Now let W denote such a closed set. A given diagram D ^ W has a total colour factor C in the 
replicated theory. This gives a contribution to the soft function in the original theory a quantity 

Y,Hd)C{d), 

D 

where we dehne C{D) to be the coefficient of the 0{N) part of C{D). This is the part of the colour 
factor in the replicated theory that contributes to the exponent of the soft function in the original 
theory. For this reason, C{D) is referred to as the exponentiated colour factor (ECF) of diagram 
D \n\- In general, C{D) will be a linear superposition of the original colour factors C{D) (this 
superposition can in principle be found by applying the replica trick algorithm as outlined in the 
two loop example). Thus, one may write 

C{D) = Y,Rdd'C{D') (50) 

D' 

for some constants Rdd'- Let us now call the set of closed diagrams W a weh, and Rdd’ elements 
of a web mixing matrix. This can be summarised as follows: 


In non-abelian theories, webs are closed sets of diagrams related by gluon 
permutations on the external lines. Each web is described by a web mixing 
matrix, which encodes how colour and kinematic degrees of freedom mix in 
the exponent of the soft function. 

It is perhaps clear that we want to consider a closed set of diagrams in some sense as a single web, 
given that these diagrams are mutually entangled. In fact, the motivation to call such sets webs 
goes further than this, as we will see in section The general structure of the soft function in 
multiparton scattering can then be written schematically as 


5 


exp 


E E J-(o)<Sc(D') 

W D,D' 


(51) 


where the first sum is over all possible webs W with mixing matrices R^^i, and the additional 
sums are over the diagrams in each web. 


It is clear, then, that the general study of webs in multiparton scattering is dominated by the study 
of web mixing matrices. These encode a potentially huge amount of physics! Namely, web mixing 
matrices encode how colour and kinematic information is entangled to all orders in perturbation 
theory, and the hope is that significant insights into the physics of non-abelian gauge theories can 
be obtained by studying web mixing matrices. 
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Figure 16: An example of a closed set of diagrams at three loop order. 


h j, k 

Replica-ordered colour factor 

Multiplicity 

i = j = k 

C(3a) 

N 

A 

II 

C{3a) 

N{N - l)/2 

i = j,k < i 

C{3c) 

N{N - l)/2 

i = k, j > i 

C{3d) 

N{N - l)/2 

i = k, j < i 

C{3a) 

nIn - l)/2 

j = k,i> j 

C{3a) 

N{N - l)/2 

j = k,i < j 

C{3h) 

N{N - l)/2 

i < j < k 

C{3b) 

N{N -1){N -2)/Q 

i < k < j 

C{3d) 

nIn -1){N -2)/Q 

j < i < k 

C{3a) 

nIn -1){N -2)/Q 

j < k < i 

C{3a) 

n\n -1){N -2)/Q 

k < i < j 

C{3d) 

nIn -1){N -2)/Q 

k < j < i 

C{3c) 

A(A- l)(A-2)/6 


Table 2: Replica trick analysis for diagram (3a), as shown in figure 16 The replica indices i, j and 
k label the gluons on the upper, right-hand and lower parts of the diagram respectively. 


4.4 A three loop example 

In order to clarify the above statements, it is useful to consider a three loop example of applying 
the replica trick algorithm, to complement the two loop example from earlier. Consider the set of 
diagrams shown in figure [l^ and let us calculate the exponentiated colour factor C of the diagram 
labelled (3a) in the figure. As in the two loop case, we must consider all possible assignments of 
replica numbers. Labelling the upper, right-hand and lower gluons with replica numbers i, j and k 
respectively, we must thus consider all possible hierarchies of i, j and k. These are given in table 
which may be compared to the two loop example of table As in the previous example, for each 
hierarchy of replica numbers, the iV-independent part of the contribution to the total colour factor 
C{3a) is found by applying the TZ operation to the diagram, and matching the result to one of the 
other diagrams in the web. One must then multiply each result by the number of distinct replica 
number assignments that are consistent with the given hierarchy. These multiplicities are given in 
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the third column of table The total colour factor in the replicated theory is then found to be 

N(N-l) 

C{3a) = NC{3a) + ^ [3C(3a) + (7(35) + C{3c) + (7(3d)] 

+ — ~ [2C(3a) + 2C(3d) + (7(35) + C(3c)]. (52) 

6 

The exponentiated colour factor in the original theory is the coefficient of the 0{N) part of this, 
which is 


C{3a) = C{3a) 


2 6 


+ (7(35) 


1 1 
~2 6 


+ (7(3c) 


1 1 
~2 6 


+ C{3d) 


1 2 
~2 6 


= - [(7(3a) - (7(35) - (7(3c) + (7(3d)]. 


(53) 


From the definition of the web mixing matrix of eq. (50), this means that the first row of the web 
mixing matrix is i(l,—1,—1,1). In fact, the full matrix turns out to be 


77 


Rdd' — X 

D 


/ 1 -1 -1 1 \ 
-2 2 2-2 
-2 2 2-2 
V 1 -1 -1 1 / 


(54) 


Problem 5. Verify this result by calculating the remaining rows of the 
matrix. 


As the above problem demonstrates, calculation of web mixing matrices in any particular case is 
straightforward, but tedious. The algorithm can be - and indeed has been - automated. Many 
more examples of web mixing matrices can be found in 77,82 , and were obtained by computer. 


4.5 General properties of web mixing matrices 

We have now seen something of the general structure of multiparton webs, and in particular the fact 
that colour and kinematic parts of soft gluon diagrams mix according to web mixing matrices. We 
have seen examples of how to calculate these matrices in particular examples. Note in particular 
that this is a purely combinatoric exercise. This suggests, broadly speaking, two different ways to 
investigate the further properties of webs: 

1. One may surmise the general structure of web mixing matrices by studying the physics of 
soft gluon graphs (e.g. the known properties of colour factors when the number of colours 
becomes large, or the need for infrared singularities to be consistent with renormalisation of 
the vertex at which all the external lines meet). 

2. One may study general properties of web mixing matrices using pure mathematics, namely 
the combinatorics underlying the replica trick algorithm. The hope is then that established 
properties of the web mixing matrices can be translated into (previously unknown) physics. 

Both of these are currently under investigation, and here we give a taste of some general features 
of web mixing matrices that have so far been established. 
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Firstly, we know that web mixing matrices have zero sum rows. That is, any mixing matrix satisfies 


Rdd' = 0, (55) 

D' 


where the sum is over any row (labelled by D). One may verify this for the two examples given 


above in eqs. (49) and (54). This result means that one may redefine the EOF of all diagrams in a 


web by a common constant. Another way of looking at this is that the contribution to the colour 
factor of a given diagram which is completely symmetric under interchanges of diagrams in the web 
does not enter the exponent of the soft function. The zero sum row rule has also been extended 
in [80] , where it is shown that certain subrows also sum to zero, corresponding to subsets of the 
web corresponding to planar and non-planar diagrams. Furthermore, a similar sum rule involving 
columns has been observed, involving extra weight factors |82j. 


Secondly, all web mixing matrices are idempotent, i.e. 

RdeRed' = Rdd'- 


(56) 


This means that web mixing matrices are projection operators, having eigenvalues 0 and 1 (each 
with some multiplicity). These properties are crucial to the cancellation of infrared singularities in 
the exponent of the soft function, as has been explored in 


55,56,82 


Both of the above properties have been proven to be completely general, including webs whose di¬ 
agrams contain three and four gluon vertices 80 . The proof of the idempotence property involves 


applying the replica trick twice in succession, and arguing that this is equivalent to applying it once 
(i.e. replicating a theory M then N times is equivalent to replicating it MN times). The proof of 
the zero sum row property uses results from enumerative combinatorics, in particular the counting 
of partitions of sets via Stirling numbers of the second kind. 


Recently, both the mathematics and physics behind web mixing matrices have been made clearer. 
In ref. [^, the combinatorics of mixing matrices was related to that of partially ordered sets, or 
posets. As the name suggests, a poset is a set of objects X endowed with an ordering operation < 
possessing the following properties: 

1. Reflexivity. If x G A then x < x. 

2. Antisymmetry. It x <y and y < x then x = y. 

3. Transitivity. If x < y and y < z, x < z. 

This ordering operation applies to pairs of elements in the set X, but such that not all pairs are 
necessarily ordered with respect to each other (hence the “partially” in “partially ordered set”). 
The relation to web diagrams is as follows: every web diagram D can be identified with a poset. 
The elements of the poset are the various irreducible subdiagrams that occur in the diagram D. 
The ordering operation < is defined as follows: given two irreducible subdiagrams x and y, x < y 
if X lies closer to the origin than y, in the sense that one cannot pull y into the origin without 
also pulling in x. As an example, consider the diagram shown in figure [T^a). In this diagram. 
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Figure 17: (a) Example web diagram; (b) Hasse diagram corresponding to the poset generated by 
the diagram in (a). 


the two gluons labelled A form an irreducible subdiagram: the fact that they are crossed means 
that one cannot shrink either gluon to the origin independently of the other one. The single gluon 
exchanges labelled B, C and D are also irreducible subdiagrams. Then the structure of the diagram 
is encoded by the information that B must be shrunk before A, C or D, and C must be shrunk 
before D. This corresponds to the Hasse diagram shown in figure [T^b): a diagram whose vertices 
represent the elements of the poset, and whose edges denote an ordering between two elements. 


Problem 6. Draw the Hasse diagrams for the web of figure 16 and hence 
show that these correspond to a (possibly kinked) chain of linked vertices. 


The language of posets and Hasse diagrams is a powerful one. Reference uses it to prove a 
number of general results, including results for the the rank of the mixing matrix for three given 
families of web, for any number of gluons. Reference carries this even further. For two particu¬ 
lar web families, solutions are given for the web mixing matrix for any number of gluons. This hints 
at the possibility that the combinatorics of web mixing matrices could be completely classified, and 
that it would be possible to go straight from web diagrams to the mixing matrix, without having 
to apply the replica trick algorithm. 


The physical content of web mixing matrices has been clarihed in ref [^. As explained above, 
mixing matrices act as projection operators, and thus have eigenvalues A G {0,1}. Associated 
with each unit eigenvalue is a combination of kinematic factors of diagrams, accompanied by a 
superposition of colour factors. These colour factor combinations have the general property that 
they look like colour factors of connected gluon subgraphs. To illustrate this, we may use the two- 


loop web of figure 14 
is 


From eq. (48), we see that the contribution to the exponent from this web 


l{C{a)-C{b)){F{a)-F{h)), 


(57) 


and the fact that there is only one combination of kinematic factors (with accompanying colour 


The colour factor appearing in eq. (57) is 


factor) is because the web mixing matrix in eq. (49) is of rank 1, i.e. has only one unit eigenvalue. 


rj-ia rjia rj-iO _ rpa rpa fjiO rj-i& _ ^ rbaCrj^a rpC rpO 

1223 1223 J 123 ^ 


( 58 ) 
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Figure 18: Connected gluon graph connecting three lines at two loops. 


where T“ represents a colour generator with adjoint index a on line i, with path ordering outwards 
from the hard interaction as indicated in figure [T^ On the right-hand side we have used the SU(3) 
colour algebra on line i: 

(59) 


rpO rpc 


= 


and we see that the surviving colour factor is directly proportional to that of figure 18, which 
is indeed a connected gluon graph. It should be stressed that it is only the colour factors that 
are related in this way: in a general gauge, the kinematic combination appearing in eq. (57) is 


completely different to the kinematic factor of the diagram in hgure [Tsl Nevertheless, this property 
is a generalisation of a similar property observed for the two-line webs of refs. 66 -68 , whose colour 
factors were found to be “maximally non-abelian”, namely equivalent to fully connected graphs. 


The fact that only connected colour factors appear in the exponent of the soft function was proven 
completely generally in ref. 81 , for all types of web at arbitrary loop order. To this end, an 


alternative formalism to the use of web mixing matrices was set up, which involved rewriting 
replica-ordered Wilson products as single exponential factors, using the Baker-Campbell-Hausdorff 
formula. Although web mixing matrices do not appear explicitly in this approach, the price one 
pays is that higher order webs must be calculated with an increasing number of effective vertices 
describing emission of multiple gluons from each Wilson line. These calculations become cum¬ 
bersome quickly as the loop order increases. However, the effective vertex formalism provides a 
well-motivated choice of basis for the connected colour factors appearing at each order 81 . Re¬ 
cently, an alternative approach to calculating soft exponents has appeared, which also involves 
effective vertices 01 0 


Interestingly, in this approach the replica trick is not used at all. 


The classification of web mixing matrices, and calculation of their associated colour factors, is only 


part of the story of understanding webs. As eq. (51) makes clear, in order for webs to be useful, one 


must also calculate the kinematic factors associated with the web Feynman diagrams themselves. 
This is the subject of the following section. 

5 Calculation of webs 


As we have seen, web diagrams are Feynman diagrams involving Wilson lines. These can be calcu¬ 
lated in either position or momentum space, and the technology needed to perform such calculations 
has developed over many years (see e.g. (^,17-19,55-64,85-^). The current state of the art is 
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Figure 19: One-loop diagram involving two Wilson lines. 


two-loops in full multiparton scattering (for both massless and massive Wilson lines), and three- 
loops for the case of two Wilson lines. 


Wilson line diagrams contain both IR and UV divergences, corresponding to taking gluon subgraphs 
to infinity, or shrinking them to the origin, respectively. One must then regulate these divergences, 
and it is conventional to use dimensional regularisation 95 in 4 — 2e dimensions. Here, however. 


a complication arises: it turns out that, for all Wilson line diagrams at arbitrary loop order, the 
relevant kinematic integrals vanish. We can see physically why this is by appealing to a one-loop 
example, namely the diagram of figure [T9l which may represent two out of many coloured particles 
coming from a given hard interaction. Labelling the momentum of line i by p*, the kinematic part 
of the diagram is (in momentum space), proportional to 



Pi 


df^k 


Pi ■ Pj 


Pj-k 


{27rY k"^ Pi ■ kpj ■ k' 


(60) 


where we have used the eikonal Feynman rule of eq. (10), and the Feynman gauge gluon propagator. 


As we saw in section such Feynman rules occur when simplifying the momentum space Feynman 
rules for whole amplitudes (i.e. with arbitrary k). In particular, we have linearised the denominators 
associated with the Wilson lines. If they were the usual denominators associated with legs of 


momentum Pi, eq. (60) would instead read 


d^k 


Pi ■ Pj 


{2-kY [{Pi + kY — m'^] [{pj — ky — m? 


(61) 


where we have taken the hard particles to be massive for full generality. One may examine the 
divergence properties of this integral using naive power counting. This essentially means that we 
identify k with some momentum scale A and see how the integrand behaves as A goes to zero (the 
IR limit) or infinity (the UV limit), ignoring the fact that all quantities are 4-vectors etc. One 
then finds that eq. (60) goes like and A*^”® in the IR and UV regimes respectively. This tells 


us that the integral is logarithmically divergent in A in 4 spacetime dimensions in the IR regime, 
which would show up as a pole in e in dimensional regularisation. On the other hand, there is no 
UV divergence, due to the terms in the denominators which are quadratic in 


One finds a different result for the linearised integral of eq. (60), which is logarithmically divergent 


(in A) in both the IR and UV regimes. Comparing the integrals, it is clear why this is: by linearising 
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the denominators associated with the Wilson lines, one removes the terms that damp the UV 
singularity in eq. (61). In other words, by taking the soft limit we have introduced a spurious UV 


singularity. A careful treatment of the integral in eq. (60) shows that these singularities precisely 
cancel each other, leaving a zero result. Furthermore, this phenomenon generalises to all Wilson 
line integrals. In other words, the UV singularities of Wilson lines are in one-to-one correspondence 
with the IR singularities of amplitudes [86] . It follows that we can extract physical information 
from Wilson line integrals by isolating their UV poles, and ignoring their IR singularities. This 
can be done using suitable regulators, over and above the use of dimensional regularisation. A 
particularly convenient regulator has been defined in refs. [^[^[^, and consists of replacing the 
usual QCD Wilson line with 


I (m) ^ 

= V exp 


i^^" / d\(3i-A{\j3i) 




(62) 


where we have shown the dependence on the coupling g, which the dimensional regularisation scale 
/i ensures is dimensionless. Also, f3i is the 4-velocity of the Wilson line, which is assumed to have 
been tilted off the lightcone if null. Equation (62) differs from the usual Wilson line due to an 


additional exponential factor, which smoothly cuts off the long-distance contribution to leave only 
UV poles in e. One recovers the conventional QCD result as m —?• 0. 


The structure of UV singularities are encoded in a quantity known as the soft anomalous dimension. 
To define this, let us first extract the UV poles of the soft function to define a renormalised soft 
function, that is UV finite: 


m 


lij,Ois{g — ) = S 




(63) 


We have here shown the dependence of each function on the regulators e and m, the coupling Og 
and dimensional regularisation scale g, and the parameters 


7*j — 


2A • 


(64) 


which are related to the hyperbolic angle, in Minkowski space, between any two Wilson lines. In 
eq. (63), Z is a multiplicative factor that collects all UV counterterms Note that this equation 
is more complicated than it looks: both S and Z are matrix-valued in the space of colour flows 


(as discussed in section 4.1), and thus the order of the factors on the right-hand side of eq. (63) is 


important. The counterterm matrix Z satisfies the renormalisation group equation 

The factor T on the right-hand side is the soft anomalous dimension, and is finite as e 
the counterterm matrix Z, it is matrix-valued in the space of colour flows. 


(65) 
0. Like 


^The fact that Wilson lines renormalise multiplicatively was first shown in refs. |88l96 
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For practical calculational purposes, there is a direct route from the unrenormalised soft function 
(in terms of webs) to the soft anomalous dimension. One may write the former as 


S ias,€) = exp 


oo oo 

E E 

_n=l k=—n 


a? 


( 66 ) 


where is the coefficient of from the sum of all web contributions at 0{a^). In terms of 

these coefficients, the soft anomalous dimension is given up to three-loop order by |82j 


r(2) = -1) _ 2 




r(3) = -6 u;(3.-i) + 




+ 3 


^( 1 , 0 ) ^^( 2 ,- 1 ) 


+ 3 


^( 2 . 0 ) ( 1 ,- 1 ) 




w 


(bO) 


^(b-l),^(b0) 


w 


( 1 ,- 1 ) 


^( 1 ,- 1 ),^( 1 , 1 ) 


(67) 


where F*^*) is the coefficient of a* in the soft anomalous dimension, and bo the one-loop coefficient of 
the QCD /3 function. We see a very general structure, namely that the coefficients of F consist of the 
single pole of the unrenormalised webs at each given order, dressed by commutators of lower-order 
web coefficients. These commutators are non-zero precisely because webs are matrix-valued in the 
space of colour flows. When this complication is absent (such as in two line scattering, or multi- 
parton scattering in the limit of a large number of colours), these commutators disappear. This is 
intimately related to the fact that webs in two-parton scatteing (in the language of refs. |66f[68] ) 
are single irreducible diagrams. Such diagrams have a single divergence at arbitrary loop order, 
such that the soft anomalous dimension (also known as the cusp anomalous dimension in this case) 
inherits all of its information from the coefficient of the single pole in each web. For multiparton 
webs, diagrams can instead be reducible, meaning that they have higher order poles in e. These 
are then cancelled by commutators of lower order webs, in a tightly constrained way. 


Problem 7. Consider the two-loop web of figure 14 Identify the single 
connected colour factor that this contributes to. Identify also the lower 
order webs that enter the subtraction term for this web, and show that they 
give rise to the same colour factor. 

The combination of an unrenormalised web with its commutator contributions has been called a 
subtracted web in refs. 55,56,82 , and has a number of special properties. Firstly, dependence on 


the regulator m used to separate IR and UV singularities cancels only at the level of subtracted 
webs. Secondly, the analytic functions of cusp angles that appear in subtracted webs are much 
simpler than for unsubtracted webs. A crucial property in making the web language useful is that 
subtraction terms operate on a web-by-web basis. That is, each closed set of web diagrams interacts 
with subtraction terms involving lower order diagrams from the same set, such that webs do not 
mix with each other under renormalisation. It is for this reason that the identification of webs as 
closed sets of diagrams is meaningful. Different webs would still mix with each other, however, 
under gauge transformations. Clever gauge choices can perhaps be used to simplify the calculation 
of webs at higher orders e.g. the conformal gauges of ref. . 
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Based on the above discussion, the procedure for calculating the soft anomalous dimension matrix 
at a given order is as follows: 

1. Choose a basis of independent connected colour factors. 

2. For each web, find (using the web mixing matrix or otherwise) the combination of kinematic 
factors accompanying each connected colour factor, if present. Note that the fact that only 
combinations are needed as opposed to individual diagrams potentially leads to simplifica¬ 
tions, due to cancellations between diagrams at the integrand level. 

3. Calculate each kinematic combination, including all relevant commutator subtraction terms. 
Further simplifications occur at this point, including independence of the regulator used to 
separate UV and IR singularities. 

4. Finally, add together all kinematic combinations for a given connected colour factor. This 
gives a gauge-invariant contribution to the soft anomalous dimension. 


This programme of work is currently being carried out at three-loop order, where webs may be 
classified according to the number of three and four gluon vertices off the Wilson lines. Diagrams 
with no multigluon vertices are known as Multiple Gluon Exchange Wehs (MGEWs). Explicit 
results have been presented for all three-loop examples connecting three or more lines 
and the two-line case has been studied in ref. 
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Certain diagrams of this type can even be 
for any number of gluons 1^. This mirrors parallel developments 


calculated to all orders i.e. 

(discussed above) in the combinatoric understanding of web mixing matrices, and suggests that it 
may be possible to understand all-order properties in the exponent of the soft amplitude. As an 
example, ref. conjectures a basis of functions which are believed to be adequate to describe 
the result of any MGEW, at arbitrary loop order. This basis, or generalisations thereof, may also 
hold for more complicated web types. An important element in the construction of this basis is 
the efficient characterisation of analytic properties using the Symbol of refs. O 


103 , which tool 


will doubtless prove to be useful at higher loop orders. Very recently, the calculation of the fully 
connected webs linking four Wilson lines has been completed 104 , in the lightlike limit. As argued 


by the authors, this is sufficient to obtain the complete soft anomalous dimension in this limit at 
three-loop order. Another interesting recent development is the use of unitarity-inspired methods 


to compute webs 57,58 


6 Conclusion 

In this article, we have provided a review of past and current developments in the classification of 
infrared singularities in scattering amplitudes. These are interesting for a variety of phenomeno¬ 
logical and formal reasons, and we have focused mainly on the case of non-abelian gauge theories, 
whose singularity structure is much richer than QED or gravity. 

Our main aim has been to provide an accessible introduction to the language of webs, namely 
Eeynman diagrams that occur in the exponent of the soft (infrared divergent) part of an amplitude, 
itself a vacuum expectation value of Wilson lines. We have seen that webs in multiparton scatter¬ 
ing are closed sets of diagrams, related by permutations of gluons on the Wilson lines. The colour 
and kinematic parts of these diagrams are entangled by web mixing matrices, whose mathematical 
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structure makes contact with interesting problems in contemporary combinatorial research. Fur¬ 
ther elucidation of these matrices, from either a physics or mathematics viewpoint, offers significant 
insights into the all-order properties of perturbation theory. This is equally true for the calculation 
of the kinematic parts of web diagrams, where systematic methods are beginning to be developed. 

Investigation of multiparton webs, and their associated physics in a variety of gauge theories, is still 
a very young subject. Many interesting results are still ripe to be discovered, and work towards 
such discoveries is ongoing! 
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